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Abstract. The single- and multi- particle inclusive spectra for strongly inho- 
mogeneous thermal boson systems are studied using the method of statistical 
operator. The thermal Wick's theorem is generalized and the analytical solu- 
tion of the problem for an boost-invariant expanding boson gas is found. The 
results demonstrate the effects of inhomogeneity for such a system: the spec- 
tra and correlations for particles with wave-lengths larger than the system's 
homogeneity lengths change essentially as compared with the results based on 
the local Bose-Einstein thermal distributions. The effects noticeable grow for 
overpopulated media, where the chemical potential associated with violation 
of chemical equilibrium is large enough. 



1. Introduction 

The theoretical study of thermalizcd hadron and quark-giuon systems is important 
for understanding of the early Universe and new phenomena in the current and 
future experiments at SPS, RHIC and LHC. The systems formed in ultra-relativistic 
nucleus-nucleus collisions produce 10'^ 10^ secondaries [0. The number of quarks 
and gluons exceeds this estimate by one order of the value if the conditions for a 
phase transition to the QCD-plasma are realized. It is reasonable to expect that 
this quasi-macroscopic system could be thermalized during small proper time tq 
after initial collision ||]. The systems formed in A+A collisions may be rather 
inhomogeneous ones even at the final decoupled stage r = r/ of their evolution 
because of a strong expansion and very high initial density. 

The current pion interferometry analysis at SPS CERN shows the effective sizes 



of such systems are i? ~ 3 7 fm 
is its geometrical mean-square size 



For a homogeneous static source the value 
, for hydrodynamically expanding systems the 



longitudinal interferometry size is approximately proportional to the hydrodynamic 
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length, i?L oc Xhydr = ^hydr l — '''/' ™ ^ Central rapidity region js). In general 
case the " pion interferometry microscope" measures the lengths of homogeneity of 
hadron systems at the final stage. Naturally, for earlier stages, t < Tf, the effective 
geometrical or hydrodynamic lengths are smaller than the mentioned ones. At 
the initional stage of thermalization the typical hydrodynamic length (longitudinal 
length of homogeneity averaged over momenta ) is Xhydr oc tq — 1 fm. 

The statistical description of a quantum-field system with small homogeneity 
regions should be done carefully. As well known the statistical hydrodynamics can 
be based on the relativistic kinetic theory Q| as well as on the method of a nonequi- 
librium statistical operator In the both methods one uses the locally equilib- 
rium distribution as the zero approximation. Then the complete Wigner function 
is usually represented by f{x) = fo{x,p) + f{x,p) and nonequilibrium statistical 
operator is p{a) = p(a) + p((j). Here fo{x,p) is a locally equilibrium distribution, 
the quasi-equilibrium statistical operator p{a) is defined on some hypersurface a 
and corresponds to maximum entropy principle under a given set of additional con- 
ditions 011^ local averages such as the energy density, the charge density, etc. The 
function f{x,p) and the operator p{a) describe the nonequilibrium flows associated 
with a heat, viscosity, etc. They give the contributions that are roughly propor- 
tional to the ratio of the correlation length (mean free path) to the hydrodynamic 
length. As to the main approximation, the distribution function fo{x,p) is cho- 
sen usually in the form of globally equilibrium distribution with the substitutions: 
T — > T{x) = 1/P{x), V Vhydr{x), etc., where the parameters depend now on point 
X. This prescription is a physically self-consistent, if the hydrodynamic length is 
much more than the Compton (or de Broglie for massless fields) wavelengths Ap of 
the quanta, Ap <^ Xhydr ■ The condition cannot be satisfied for the effective mass 
TO < 0.1 GeV at the early stage, t ~ Xhydr ~ ^f'm , of the matter evolution in 
nucleus-nucleus collisions. One of the aims of this paper is to find the locally equi- 
librium function fo{x,p) in the general case of an arbitrary ratio Ap to Xhydr- We 
shall use, for the purpose, the method of the quasi-equilibrium statistical operator 

0, @, i- 

The another problem, we discuss here, is how to calculate the two (many)- parti- 
cles inclusive spectra for very inhomogeneous locally equilibrium systems. Strongly 
speaking, it is impossible to define the one-particle Wigner function in this case: 
the necessary condition that the thermal averages < {p+ ^)a{p — ^) > have to be 
diagonal enough, i.e. q^yj = A^^^^ ^ TO^[|j, does not satisfied. In such a situation 
the direct use of the method of statistical operator is appropriate for the spectra 
calculation. 

The special interest for discussion is the appearance of the additional terms 
in the double particle inclusive spectra connected with the non-zero averages of 
the products of the creation and annihilation operators, (a+a"'") and (aa). Such 
terms have been obtained first for pious in the Gaussian current model (l4|. The 
additional terms in the correlation function do not appear in the nonrelativistic 
quantum- mechanical approach m^. We shall consider all these problems using the 
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generalized Wick's theorem for thermal locally equilibrium systems and give the 
explicit analytical structure of the pair-correlation function for different sorts of 
bosons. 

In the section 2 we deal with the description of inclusive spectra in the standard 
Wigner's representation and discuss some basic points of this approach. 

In the section 3 we develop the statistical operator formalism for calculation of 
bosonic operator averages, (a+(pi)a(p2)), (a(pi)a(P2)) , etc., in the locally equilib- 
rium systems and discuss the physical conditions for the hydrodynamic solutions 
we are interested in. 

In the section 4 the thermal Wick's theorem is generalized for locally equilibrium 
systems. This is the basis for the calculation of the double- and multi- particle 
inclusive spectra. 

The section 5 is devoted to the analytical calculation of the single particle spec- 
tra for a physically important case of the boost-invariant expansion of a hadron 
and/or quark-gluon matter. We derive there the correction term to the Bose- 
Einstein heat spectrum and demonstrate the tie of the term with the spectrum 
of the so-called Milne's particles. 

In the section 6 we obtain the structure of the pion-, kaon-, and photon- pair 
correlation functions. We calculate also the maximum value of the interferome- 
try peak for these particles and the analytical approximation for the two-particle 
correlation function in typical experimental situations when mr/ ^ 1. 

2. The Statement of the Problem 

The description of the inclusive spectra and correlations for a multiparticle produc- 
tion is based on a computation of the following type of the averages 

where and ap are the creation and annihilation operators, corresponding to a 
quantum field ipout{x,t) when an interaction is switched off. The brackets <...> 
mean the average over some density matrix describing the state of the system on 
a some hypersurface a. In the S'-matrix theory the state is the |oiti)-state at a : 
t = oo. In the statistical thermodynamic models of a multiparticle production the 
density matrix is chosen to be the statistical operator p and surface a is usually a 
freeze out hypersurface. The averages (|l|) taken on this hypersurface are coincided 
approximately with ones taken on an arbitrary hypersurface that situated within of 
the light cone of the future as to the freeze-out hypersurface. It corresponds to the 
preservation of the momentum distributions of free streaming particles in Eqs.([l]) if 
one neglects the final state interaction and Coulomb corrections. The hypersurface 
a can correspond also to an earlier stage of the evolution if one studies the dilepton 
or photon productions from a hadron and/or quark-gluon plasma. It the latter case 
the operators and ap have to correspond to weekly interacting quasi-particles 
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with the standard relativistic form of the dispersion relations in the medium and 
the masses depend now on the temperature and density. In this paper we shaU not 
consider such a situation in details. 

The inclusive double particle spectrum in (Q) is usually calculated under sup- 
position that the four-operator averages can be decomposed into the products of 
the irreducible two-operator ones 

(a+ a+ Opiflp,) = (a+ flp,) (a+ a^,) + (a+^ap^) (a^^ap,) (2) 

In this case the problem of the inclusive multi-particle spectra is reduced to the 
calculation of the (ap^flpa) averages. 

If one considers the free identical particles 1 and 2 the two-operator average 
can be expressed by means of the following distribution function 

f{x,p) = (27r)'3 J ^4^5(p . „ ^) ^ 

where p — {pi + p2)/2 and the average is done in a space-time region where the 
interaction is negligible. Indeed, if we consider some hypersurface S that situated 
in this region and can be closed by a plane surface t = const, it is possible to use 
the equation 

J dY^^p'^e'^-"' = (27r) Ve'''°*<5^(k) if p • fc = (4) 
that follows from the Gauss theorem. Then 

{a+{pi)a{p2)) = J d^up^e''^'i"-^"'^H{p ■ u)5^{ci - u) {a+{p ~ |) a{p + §)) 

(5) 

Here a is the part of the hypersurface ^ where f{x,p) ^ 0, q — P2 — Pi- 

The expression (|^) is general and describes the operator's averages for the radi- 
ating matter when the hypersurface cr is an arbitrary hypersurface situated within 
of the light cone of the future as to the decoupling 4- volume. In the general case the 
function f{x,p) is rather complicated, even not positively defined. The distribution 
function f{x,p) is coincided with the single- particle Wigner function fw{x,p) for 
free fields if < a^{p + |)a(p — |) > is diagonal enough: q^^j ^ m^. Then there is 
the direct tie between the distribution function (|^) and complete Wigner function 
N{x,p) = i47r(27r)-^/d4wexp(-ipu) <: Lp{x + ^v)Lp{x - ^v) :> p. 

The real calculations of the final spectra and correlations simplify greatly if a 
system is thermal and decoupling volume is narrow enough in time- like direction and 
so it can be considered as the freeze-out hypersurface. In this case one can use the 
thermal matrix density p at this freeze-out hypersurface a and calculate the phase- 
space distribution function f(x,p) directly. When the surface a changes from event 
to event it is necessary to do the additional average over a in the all final expressions 
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for spectra. But this procedure cannot be used even formally for description of two 
( many)-particle spectra when the radiation volume is an essentially 4-dimensional 
one. For this aim instead of single particle Wigner function (let us suppose here that 
/(x,p) = fw{x,p)) one have to use the quasi-classical density of particle emission 
g{x,p) = p°(fN/d^xd^p. The latter can be expressed by means of the derivation 
of complete Wigner function, p'^'-df^N{x,p), that take into account the interaction 
in the system leading to continuous radiation during some finite time. The decay 
of the resonances is one of an example of such a 4-volume emission. We will not 
consider this case in the paper. 

As we show hereinafter, it is convenient to choose the one-component scalar 
field as the basic model of our consideration. The Lagrange function has the form 
L=Ijkg+Lint corresponding to the free Klein-Gordon field and an interaction 
term. If the latter is characterized by the coupling constant a (with the dimension 
equal to product of energy and volume), we can neglect the interaction energy and 
momentum if the following conditions for temperature T and particle density n{x) 
are satisfied Q 

^ « 1 (6) 

Neglecting the interaction terms, the single-particle Wigner function near the 
mass-shell is associated with the local current To be simple we will consider 
here the real field and use the particle fiow as the current 

j^(x)^^^+\x)^^^'H^) (7) 
where the decomposition of the field into "positive" and "negative" parts looks like 

^(x) = v^+\x) + ^^-\x) EE [2{2n f]-i J ^ [a+e^^-^ + a^e'^^'^] (8) 

The expressions for the single particle spectrum and the double particle corre- 
lation function follow immediately from Eqs. (|l|) , (||) , (|^) : 

= j da^p^f{x,p) (9) 



C{pi,p2) = 1+ [pIpI 



da,,p^^e'^-^f{x,p) 



(10) 



If one considers the infinite homogeneous thermodynamic system (with the 
chemical potential /i, the energy-momentum operator P and the operator of par- 
ticle number JSl) that moves as a single whole with 4- velocity the result of the 
averaging over the equilibrium statistical operator 



1 

Peq = exp 



P^u^ + Hm/T (11) 
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is the following |^ 

(a+(pi)a(p2)) = (2^)V/s.iJ.(p)^'(Pi - P2), {a{pi)a{p2)) - (12) 

where Jb.eXp) the Bose-Einstein distribution for the globally-equilibrium sys- 
tems. Let us put for simplicity fi = 0. Then 

fip,x) = 7^-^^ 7 = fB.E.ip;P,u); (3^1- ^ const, u" = const. (13) 

exp(pp • u) — 1 J 

Here the Wigner function does not depend on x. The thermal Wick's theorem 
can be proved for a such type of systems, that leads to the result (||). The main ap- 
proximation for the Wigner function of an expanding hadron and quark-gluon gas is 
usually based on the distributions like ( [l^ ) with the substitutions (3 = const — > /3(a;), 
u = const — > u(x). Such substitutions are physically reasonable if the wavelength 
of the corresponding quanta, Ap — l/rne//, is much less than the hydrodynamic 
length 

Ap << Ahydr oc min{|9^u* I ^ ,\gradT\^^} (14) 

Asterisk marks the values in the (local) rest system. As it was mentioned in 
Sec.l, the typical hydrodynamic length is approximately equal to the proper time of 
the hydrodynamic expansion, Tf , and the inequality (|lj) is satisfied if mr 3> 1. For 
pions and kaons mrj >1, for chiral quarks and gluons rnr ^ 1. For thermalized 
photons there is always the momentum region where inequality ( p^ ) is strongly 
violated. 

At the end of this section we would like to emphasize that the results ([T^), 
( p^ ) for averages (a+ Op^) as well as the zero values for the averages {a-p-^cip^) have 
been derived for infinite homogeneous equilibrium systems only and cannot be auto- 
matically applied to locally equilibrium inhomogeneous system by using the simple 
substitution P = const P{x), u = const u{x) in the Bose-Einstein distribution 
(p^. This concerns also of the two-particle spectra ( p^ ) and expansion (|^) which 
is based on the Wick's theorem for globally equilibrium systems. 

It is interesting to mention that despite of the different structure of the correla- 
tion function (p^ ) in different approaches (e.g.,[^o|-0) they will give approximately 
same results being applied to systems that are quasi- homogeneous ones, mrj 3> 1, 
and are described by the same Wigner functions. But for strongly inhomogeneous 
systems it is impossible to introduce by a standard way the single particle Wigner 
function as well as to preserve the structure of the correlation function (|l^) based 
on Eq.(^. Therefore, all these approaches are off the region of their applicability. 
One have to calculate the averages such as {cLp^Up^'^, {^a^_^ap^'j directly. The formal 
distribution functions defined by (^) will be differ from local Bose-Einstein distribu- 



tion even for ideal Bose gas. The structure (10) of the correlation fimction will be 
destroyed altogether with the Bose-Einstein distribution (^3|). In the following sec- 
tions we propose the method for study of spectra and correlations in inhomogeneous 
thermal systems. 
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3. The Method of Locally Equilibrium Statistical Operator 

The hydrodynamic description of quantum-field system, as known, can be based on 
the method of non-equiUbrium statistical operator ||^ , Q , |^ , |l6j . The initial step 
in this method is to construct the so-called quasi-equilibrium statistical operator 
that describes hydrodynamics of the system neglecting the viscosity effects, heat 
conductivity, etc. In other words, the operator describes the locally equilibrium 
systems and since it will be used for this aim only we will call it as the locally 
equilibrium operator p. To build the operator one usually applies the maximum 
entropy principle 0, [|l^. The method carries into effect in a full analogy with 
the Gibbs method for homogenous equilibrium systems. In the latter case the set 
of the averages < E >, < P >, < Q >, etc., is considered as the fixed additional 
conditions when the entropy S is maximized. For the locally equilibrium systems 
considered on some hypersurface a with a time-like normal vector n'^ the collection 
of the additional conditions is based on densities of energy £(a;), momentum p(a;), 
charge q{x), etc. In the relativistic covariant form they look like 

(n,(x)f^'^(x)) , in,{x)J-{x)) (15) 

where T'^'^{x) is the operator of the energy-momentum tensor, J'^{x) is current 
density operator. According to the definition 

p =e-'5('^) (16) 

The entropy is maximized under the additional conditions like ( p^ by the 
Lagrange multipliers method 



S — max Sp 



—p\np—p / da n^{x)[PyT'^^ {x) — pf3J'^ (x)) — Xp 



(17) 



where Lagrange multiplier P^{x)— u'^ (x) /T{x) ||16|. The formal solution of Eq.(|l 
gives us the result for entropy M , flffl 



S = S{a) = $((t) 



da n^{x){(3^T'"'{x) - pP.r{x)) 



(18) 



where $(cr) = In 5pexp{J da n^{x){f3^T'^^{x) — p[3J'^{x))} is Masier-Planck func- 
tional. Let us put p = Q for simplicity. Sticking to the analogy with the method of 
statistical operator for globally equilibrium systems where all the operators under 
the sign Sp are mutually commuted, we demand that 



n,{x)P^{x)f^^''{x) , n^{y)Ps{v)f'^{y) 



= 



(19) 



In this case the set of the additional conditions ( [l5| ) has the standard inter- 
pretation and we hope to avoid some mathematical difficulties that may possibly 
appear in a more general case. 
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As we mentioned in Sec.2, we shall start from the free one-component scalar 
field (||) with the standard commutation relations for the operators 

[a{p),a+{p')]=p''dHp-p') (20) 
The energy-momentum tensor has the form 

The commutation equation ( p^ ) is solved using the Eqs.(^, (^), ( pi] ) and gives 
the following conditions for hydrodynamic values taken on a hypersurface a : 

(22) 

X ■ u{x) — y ■ u{y) \fx, y <E a 

It immediately follows from the Eqs.(^2|) that the commutation relation (|l^) is 
satisfied when one of the following conditions is realized: 

1. The hypersurface a is plane and u^{x) = const {v{x) = in the reference 
system where a is t = const ) . Actually, this means that the system occupies the 
space-time region with some distribution in temperature (falling down from central 
highly excited part to the vacuum at the periphery) and has no internal motion. 

2. The hypersurface a is defined by the condition — x| = = const, 
and the one-dimensional expansion along L-axis occurs with 4-velocity uq = t/r, 
ul = xl/t, ut = 0. If /3 is a constant on a, it is reduced to the well-known 
boost-invariant expansion , , that is the basic model for application of the 
hydrodynamic theory to multiple processes in high energy collisions. 

3. The hypersurface a is defined by the condition t^ — x^ — t"^ — const and 
there holds the two-dimensional expansion uq = t/r, \it = xt/t, ul — 0. 

4. The hypersurface a is defined by the condition — x^ = = const, 
3-dimensional hydrodynamical expansion has the form u^^ — x^ jr. 

The method proposed to be used to find the averages of the operator products is 
based on the Gaudin's idea [l9| ] for globally equilibrium systems and is the following. 
We represent the locally-equilibrium statistical operator p defined by (p^), ( p^ ) in 
the form 

^ ^ da,(3,,f''''{x) 



P=^exp 



(23) 



where the integral is taken over corresponding hypersurface a as it was discussed 
before. Let us introduce the operators that dependent on some parameter a in the 
following way 

(24) 
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and use the matrix notation 

It is easy to get the fohowing equations 

A(p,a-0)-A(p) (26) 
< A{p)a{j)') >=< a{p')A{p,a = 1) > 

(27) 

< A{p)a+{p') >=< a+(p')A(p, a = 1) > 

The latter equations follow from the trace invariance under the cyclic permu- 
tation of operators. To express the operators A{p,a) through a+(p) and a{p) we 
shall use the equations that stem directly from Eq.(E4) 



da 



Mp) , / da^P^T'^'^ix) 



-a J da,,l3,,T''''(x) ^^^^^ 



Using the commutator (|2y) and conditions (^2|) we find the concrete form of 
Eq.(|8l) for scalar field (|) 



dA{p, a) 
da 



(fkK{p, fc)A(fc,a) 



(29) 



where the matrix kernel of integro-differential equation ( |29| ) has the form (asterisk 
means the complex conjugation) 



K(p, k) 



G{p,k) G{p,k) 
-G*{p,k) -G*{p,k) 



(30) 



where 



G{p,k) = -j^Jda ^k-pyx^p. [k-up-u-{k-p--m'^)/2\ 
G{p,k) = jda e-^ik+pyxHp, [k . u p ■ u - (k ■ p + m^)/2] 



(31) 



Here f3{x) is the inverse of the local temperature T{x), u{x) = n{x) is the hydrody- 
namic 4- velocity and we use here the integral measure in the form da^ — da n^{x). 
The solution of the system of integro-differential equations (p9|) , which is con- 



sidered according to Eq.(26) as the Cauchy problem at A{p, a = 0) = A(p), is (see 
Ref.g): 



A(p, a) = A{p) + Y.^J d'fcK„(p, k)A{k) 



(32) 
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where K„(p, k) is the n-th iteration of matrix kernel K: 



(33) 



K„(p, fc) = / dsidS2—dSn-l'K{p, Si)K{si, S2)...K{Sn-l, k) 

Using Eq.(^^ we have the properties for the iterations of the kernel 

Kf = (-1)"K:", = (~1)"K:12 (34) 

The system of integral equations for the operator averages follows from the 
solution of the integro-differential equation ( p9| ) and the relation for averages 

Er=i I d'k {Ki\p, k) {a{p')a+{k)) + K^p, k) {a{p')a{k))) = -p^d^P p') 
En=i I d'k {Kfip, k) {a(j>')a+{k)) + Kf{j>, k) {a{p')a{k))) = 

(35) 

and 



T.n=i h I d'k {KHp, k) {a+{p')a+ik)) + K^{p, k) {a+{p')a{k))) ^ 

EZi i I d'k (K2i(p, k) {a+{p')a+{k)) + Kf{p, k) {a+{p')a{k))) ^ p^S^p ^ p') 

. _ _ . . . (36) 

The integral equations ( |35|) and ( |36| ) contain the complete information about 
the one- and many-particle inclusive spectra for the locally equilibrium thermalized 
Klein-Gordon field. If the system is an infinite homogeneous (/3=const ) one and is 
considered on a flat hypersurface a: = const in the rest frame where = const, 
it immediately follows from (P) (^ that K" = -Pp°S^{p* - k*), (po ^ p ■ u), 
and according to Eqs. (35) and (^) we have 



(a+(p)a+(p)) =(a{p)aip'))^^^0,{a+ip)a{p')) = ^"fJP ^\ (37) 



eg \ /eg 



exp(/3p • w) — 1 

This corresponds to the standard result (|lj) for globally equilibrium systems 
and leads to the Bose-Einstein distribution for a homogeneous ideal gas. In all other 
cases the result will be different from (p7|); the averages (a+a"*") and (aa) do not 
vanish because c-factors attached to the corresponding operator pairs in the tensor 
T^^ do not become zero after the integration over a. Hereinafter we shall consider 
the solution of Eqs.(p5|) and (pfl) for the concrete locally equilibrium systems. 



4. The Thermal Wick's Theorem for Locally Equilibrium Sys- 
tems 

The double and multi-particles inclusive spectra are defined by Eqs.(^ and are 
expressed through the four and many operator averages. If a system is in globally 
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equilibrium state, there can be used the thermal Wick's theorem to express the 
many-point operator averages as the products of two point ones. For 4-point average 
the corresponding expansion is given by the Eq.(^. Our task now is to generalize 
thermal Wick's theorem for locally equilibrium systems. 

Let us consider the 4-point operator averages. First we introduce the notation 
for the operator expression 



oo 

E 

Tl=l 



^ j d?kKn{p,k)^K{p,k) (38) 



in order to simplify computation. Then the integral equations (^5|) and (pq ) take 
the compact form 

-/(53(p-p') 



K{p,k) (a(p')A(fc)) = 



K{p,k) (a+(p')A(fc)) 






V°5\V - P') 



(39) 



Using the trace invariance under the cyclic permutation and solution ( p2| ) one 
can get convinced that 



(A(pi)a+(p2)a(K)a(p'2)) = {a+{p2)a{p\)a{p'^)A{pi,a ^ 1)) = 
{a+{p2)a{p'^)a{p'^)A{p{)) + k{puk) {a+ {p2)a{p[)a{p',)A{k)) 



(40) 



After the commutation of the vector A and the representation of the arising S- 
function by means of the Eq.(^9|) we have 

(A(pi)a+(p2)a(pi)a(p^)) - {a+{p2)a{p[)a{p'^)A{p,)) = 

K{pi,k)[{a{p[)A{k)) {a+{p2)a{p',)) + {a+{p2)A{k)) {a{p[)a{p'^)) + (41) 

(a(p^)A(fc)) {a+{p2)a{p[))] = K{pi,k) {a+{p2)a{p[)a{p'^)A{k)) 

Just in a similar way one can derive the analogous equation using the substi- 
tution a'^{p2) a.{P2) and the Hermitian conjugated to it. It is also worthy to 



mention that the equality between the last two parts of Eq.(41) is preserved when 
the vector A commutes to the left side in all the brackets. The Eq.(^l]) is easy gen- 
eralized for the case of any even number of the operators by the induction method. 
All that means that the following equation is valid for even number of operators 

k{p,, k)A {A{k),A{p2), ...,A{pj), ...) = (42) 

where 

A = {A{k)A{p2)...A{p,)...)-Y,\{A{k)A{p,,)) n {A{p,A{p^,))\ (43) 
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and A{p) is a+(p) or a{p) and V is the permutation sign. If the integral operator 
if is a nondegenerate one, the Eq. (^3|) has the unique solution A = that means 
the average of any even number of operators expands in the sum of the products of 
all operator pairs taken in the same order as they were in the initial expression: 

{A{p,)AiP2)-A{pn)) -EH (^fe^fe')) ' (" = 2fc) (44) 

V ]'>] 

It is obvious that the averages of odd numbers of the operators are equal to zero 
because of the bilinearity of the energy- momentum tensor T^'^{x) in operators a"*", a. 
So, the theorem is proved. The main peculiarity as compared with the standard 
results is the presence of additional terms like (a+(pi)a^(p2)) and {a{pi)a(p2)) in 
the expansion. As it will be shown in the next section these non-zero terms arise 
because of a space-time finiteness of the homogeneity regions in locally equilibrium 
systems. 



5. Boson Spectra in the Boost-Invariant Hydrodynamic Model 

In this section we consider the non-trivial case of a locally equilibrium system sat- 
isfying the conditions (|2^). It corresponds to the well-known and the physically 
important hydrodynamical solution of the ID boost-invariant expansion p^ . 

Let us introduce the standard variables to analyze this hydrodynamic solution. 
The space-time variables in terms of rapidity y and proper time of the expansion r 
look as follows 

t = T coshy, XL = T sinhy, = coshy, ul = sinhy, dcr^ — u^TcPxTdy (45) 

The system is considered on the hypersurface t = const, where the inverse of 
temperature (3 = const. 

The particle momentum can be also expressed in terms of the particle longitu- 
dinal rapidity: 

p = (mp cosh 6'p, pt, TOp sinhSp) , wip = tot [p) — \/ 'm? -\- 

(Pp = nip coshOpCpprdOp] (46) 

[a{p),a+{p')] = [a(pT, 0p),a+{p'j., Bp,)] - ^^(p^ _ p^)^(0^ _ g^,) 

To simplify the problem let us assume the transverse radius of the hydrodynamic 
tube to be much larger than the hydrodynamic length t. So one can neglect the 
influence of a finite transverse size of the system on the form of spectra. The 
calculation of the basic functions (pT|), G and G, that is easy to do using the 
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variables (Bq) with an infinite transverse region, give us the following results 



oo 

27r cosh 6k J 



i2m}~TZ sinh ( — 



2 



G{p, k) = ^I^l^S^ipr + kT) [ dzVz2-lexp 
TT cosh Ok J 



'i2mkTZ cosh 



(47) 

Ok ' 



_ _ (48) 
The functions (|47|), (|4q ) are the distributions (the generalized functions). For 
example, for globally equilibrium systems G(p, k) cx (5'^(p — k). They can be con- 
sidered as the Fourier transforms of the distributions vT+z^ and 9{z — z'^ — 1 
acting in the space of rapidly decreasing functions f{9). According to the general 
rules of the operations with the so-called tempered functions , we shall mean or 
directly substitute the regular functions 



Ge = G ( VT+^ e-^l^l VT + ^) , G, = G ( - 1 ^ e^^l^l \/z^ - 1 ) (49) 



instead of {]^, ( |4§| ) and use limit e ^ in the final expressions. 

Let us rewrite the integral equations ( p5| ) in the variables (^), (p6|). Note that 
the element K*f of the n-th iteration of the matrix kernel consists of even number 
K^^ or K^^ in Eq. (js^) if i = j and odd number of them if i ^ j. So taking 
into account the structure of the operators G and G, the solution of the integral 



equations (35) can be presented in the form 
{a{p)a+{k)) ^ Ai{9p-9k)6\pT--kT), {a{p)a{k)) ^ A2{9p- Ok) S^pr + kr) (50) 



Then we have after the integration of Eq.(|35|) over tPkx 

f^-J d9kKr,{9p - 9k)A{9k - 9^,) = ( ^^'^ ) (51) 

n=l ^ ^ 

here the matrix A = {Ai,A2) is defined by Eq.(pO[) and K„(0p — 9k) are the n-th 
iteration of the kernel 

'^^''-''^-[-G*{ep-9k) -G*{9p-9k) ) ^''^ 



The matrix elements in ( p2| ) are defined by Eqs. 

G{9p - 9k)5^{pT - kT) = mT cosh ^^fcG(p, k) 

(53) 

G{9p - 9k)5''{pT + kr) = mTC0sh9kG{p, k) 
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Fourier transform of Eq. (|51|) gives 

[expK(t)-I]^(t)= ) (54) 

where I is unit matrix, K(i) is Fourier-transformed matrix js^). 

The solution of the Eq.(p4|) is easy to find after the diagonahzation of the K(t)- 
matrix. The proper values of the corresponding characteristic equations are 



-u;{t), UJ2 = Lo{t) = \]G^{t)-G*{t)G{t) (55) 

and the matrix U transforming the coordinates Ai and A2 of the vector A to the 
new ones A!^ and A!2 and with the diagonal matrix K has the form 

Now one can get the solution of the Eq.(p^). It is the following 

/ + \ A ^ 1 (G - w)2 + G*Gexpa; 

{a+a)=Ai-l = (57) 

^ expcj-1 [G-ujf-GG 

I X A G* (G - Lj) l + coshw 

aaV =A2 = ^ =j= (58) 

* (G-w)2-GG smhcj ^ ' 

Before we will analyze the analytical structure of the Fourier components (p7|) 



and (58) of the pair operator averages, it is necessary to consider the principal 
problem. According to the construction of the matrix density p ( p3| ) the physical 
vacuum defined as a |0) = is not the proper vector of p because of the presence 
of the terms proportional to a+a+ in the energy-momentum tensor. The last term 
does not vanish after the integration over da except for the case a: t — const, 
/3 = const. This means that at zero temperature, in the limit /3 ^ 00 (T — > 0), the 
operator averages such as (|5^ ) and (^) do not describe the modes of the physical 
vacuum |0) but is associated with the modes of some "lowest" state |0') of the 
operator J da^(3^T^^'^ {x). So all the exciting modes at finite temperature appear 
over background state |0') and we have to renormalize the averages: 

(a+a);;™ = {a+a)^ - {a+af^^^ - (cxpcc; - l)-i(l + 2 (a+a)f^°°), etc, .. (59) 

After this important remark we have finally 



/ _i \ren 1 |G| ren 1 G , i ,^ren 1 G 

(a a) — , (aa/f = , (a a ) = 

' exp uj — 1 uj exp uj — 1 uj ' exp uj — 1 lo 

(60) 

where uj is defined by Eq.(p5|) and 
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G{t)^- / dxVTT^ J d0e»2rn^r . sinh 0+ ^2te 



9 OO 

G{t) = lim,_o / d0e'*^ J dzV^2^g-<2m^ 



r 2 cosn TT—e z 



(61) 



(62) 



= imrP J d9,^^^Hf'{2mTT cosh 



(2) 

Here K^, is the modified Bessel function and H[ '{z) is tire Hankel function 
of second order. The functions G{t) and G{t) exhibit the following asymptotic 
behavior: 



rriTT ^ 1 



G{t) - -(3mTWl + tyimrrr + ^^^^ (l + t^/{mTTf) 



\G[t) \ = f3mT 



1 



2mTT 



{l + tyiniTTf) ' + 0{{mTT)-^) 



uj(t) ^ PrriT l + t^/{rnTT) 



mj'T <C 1, i < Tot''" 



2 1 - 3t/2mTT - t^/2mTT 



6(TOTT)2(l+t2/(TOTT)2)' 



(63) 

(64) 
(65) 



G{t)^ 



1 



, oITLtT , TOrT 

In^— !— +ln— ^ + 1 



+ O ((tott)^ 



G(i)---^{l+(^[ln^n^+lnn^ + l-.r2]} + 
In 12X1 + o ((tott)4) 



(66) 



(67) 



(68) 



The asymptotics is found by the saddle-point method for to^t ^ 1 and using 
the representation of Eqs.(|6l|), (62) as the sum of the hypergeometrical functions 
(see [Q) for tott < 1. 

It is important to note that function lo is described by the expression (^8|) for 
all values TOtt with the accuracy more than, at least, one percent as the numerical 
calculations have shown. It means, that the main factor (expo; — 1) ^ in Eqs.dl) 
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is actually the Bose-Einstein function (^3|). We believe that there is deep analytical 
foundation of this fact but wc cannot explain this mathematical gift now. 

It is worth to mention that the introduction of a non-zero chemical potential 
associated with violation of the chemical equilibrium in boson gas (overpopulation) 
leads to changing of all results by the substitution: G{t) G{t) — /i. 

Let us calculate the single particle spectrum for the case /i = 0. The pair 
operator average is described by the inverse Fourier transform of Eg. (|60|) . With 
Eq.(^0|) taken into account we have 

{a+{p)a{p')) = 5\pT - p't)^ J (a+a)^ (69) 

(a(p)a(p')) = {a+{p)a+{p')y = S^pr + p^)^ J dte-^^"^-'^'^' {aa), (70) 

One can calculate the formal distribution function defined by (^) using the 
variable C = Arsinh 2mr"cosh e ^^'^ representation ( ^6| ) for momenta. 

(^a~^ aj (t—niTT sinh(6' — y)) 



(71) 



(2^ 

The hydrodynamics rapidity y defines the position of point x according to 
Eq. (^5|) . The last approximation in Eq.(^l|) is valid when the Fourier transform 
F [{a~^a)^] (C) is fast decreasing in (, that is satisfied. So in the local rest system of 
the fluid element, y — 0, the result has the form 



rriTT ^ 1 



d^|u=o " ^^'^'^^'^1"=° exp(/3po) - 1 \ + ^ ^ 



ttitt <C 1 



-/(X,p).|u^0« , ^fBE{p){2fM{p)) (73) 



d^xd^P\^^Q ' exp(/3po) - 1 Trrpo 

We introduce here the designation 

/m(p) = (a+a)f-*°° ■ . « + J = (e'""^" - (74) 

The function /m(p) describes the " heat" spectra with the "temperature" Tg// = 
l/27rT and the real temperature T = for the so-called Milne's particles ||2^], which 
appear when the 2-dimensional field system is quantized in the hyperbolic space- 
time that is known as Milne's Universe [p4[P]. It happens due to the mixing of the 

^The same spectrum with t = 1/a is well known also for the Rindler's particles that appear in 
uniformly accelerating reference system with acceleration a due to the quantization on a time-like 
hyperboloid formed by the world lines of accelerating observers. 
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positive and negative frequency components of a field in the hyperbolic world in 
compare with the Minkovski one. In the boost-invariant hydrodynamic model this 
hyperbolic space-time is formed by the isotherms. The state of the "lowest" energy 
for the operator J doj^u^T'^'^ containing in the statistical operator p was found in 
the Ref. in the limit corresponding to r — > in our case. So the unrenormalized 
spectrum {a^a)^^°° gives the standard result (|7^) for the spectrum of the Milne's 
particles in this limit. After the background subtraction procedure, which is neces- 
sary since we study the particles against the background of the Minkovski vacuum 
but not the Milne's one, the trace of this phenomenon can be observed at the fi- 
nite temperature only as it is demonstrated by the Eqs.(^), (^). The physical 
reason for this lies in the creation of additional quanta of a bosonic field due to the 
interference of the positive and negative frequency components of a field when the 
latter begins to embrace the " hyperbolically" expanding medium. Naturally this 
effect may be noticeable only if the wave-length of the quanta is much large than 
the length of homogeneity , in other words, when non- localized quanta "feels" a 
space-time inhomogeneity of a medium. 

The numerical demonstration of this effect for the distribution function is rep- 
resented in Fig.l. 




Note, that at ^ there is noticeable deviation of the distribution from the 
Bose-Einstein result even at rriTT 2> 1 that began to be very large when /x n^r = 
m(l — l/2TOr) : 

fBE{p,x) { (2rpo(Po-A^))V ^ ' 

The single particle spectra according to (|l|),(|69|) for large enough radius R of the 
hydrodynamic tube has the form 
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/ 



(76) 



The effects of inhomogeneity will be important at the early stage, 1 fm for 
low-mass gluons or quasi-bosons if such objects there are at this stage. It could 
lead to the enhancement of photons and dileptons with low invariant mass. At 
the final freeze-out stage, even if uitt ^ 1, the distortion of spectra due to the 
inhomoheneity will take place when the violation of the chemical equilibrium in 
boson gas is strong: the chemical potential /i is closed to /icr- 

Note, however, that we did not take into account a longitudional geometrical 
size of the system, R^.H means that our consideration is limited by the condition 
p° > 1/Rl- The condition p° ^ 1/Rt have been supposed before. Note, that as it 
was found in a smallness of transverse geometrical size of the system, Rt < 
leads to an reduction of number of soft quanta in locally equilibrium system in 
comparision with the Bose-Einstein spectra at the same values of thermodynamic 
parameters. 

6. Two-Particle Correlations 

The generalized Wick's theorem (^4|) allows one to consider many-particle inclusive 
spectra and correlations for concrete fields. Let us begin from pion interferometry. 
The isotopic invariant Lagrangian of free pion field in the representation of real 
3-component pion field operator vr and nonzero commutators have the form {i,j = 



So it is easy to see that dynamic equations split into three independent ones for 
each real component ttj ; the statistical operator p is the product of three commuting 
exponents and all the previous results are preserved for each Hermitian field tTj 
ip independently. 

As it well known, the creation and annihilation operators of pions tt^ , 7r_ , ttq are 
described by the pion field in the complex representation and are connected with 
corresponding operators of the Hermitian (tti , 7r2 , tts )-field in the following manner 



1,2,3.): 



,x/) =-^;.^(x-x' 



(77) 




The results are 
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• for vr vr (and similarly 7r+7r+ ) pion pairs 

{aX{pi)aX{p2)a+{pi)a+{p2)) = {a+{pi)a{pi)) {a+{p2)a{p2)) 

+ {a+ {pi)a{p2)) {a+ {p2)a{pi)) 

• for vr+TT^ pairs 

(aJ(pi)at(p2)a+(pi)a-(P2)) = {a+{pi)a{pi)) {a+ (jp2)a{p2)) 

+ (a+(pi)a+(p2)> (a(p2)a(pi)) 

• for vr^TT^-pairs 

{at {pi)at {V2)m{pi)m{p2)) = {a+{pi)a{pi)) {a+{p2)a{p2)) 
+ {a+ipi)a{p2)) {a+{p2)a{pi)) + (a+(pi)a+(p2)) (a(p2)a(pi)) 



(80) 



(81) 



(82) 



Such a structure of the two pion spectra has been gotten in the Ref.|T^] in 
the quantum optics model (Gaussian random sources) and in the Ref. p5|] in two 
particle quantum mechanical approach. The significance of our result is that we 
calculate the averages (a'^a) , (a"'"a^) and {a a) in the concrete realistic model of 
locally equilibrium hadron medium. The correlation function C{pi,p2) for different 
pion pairs can be obtained by the division of the two particle spectra (|80|)-(|8^) by 
the single particle momentum distributions (^9|) of the corresponding pions. The 
theoretical maximum of the correlation functions C{pi,p2) is achieved at Pi = P2 = 
and can be easily estimated using the saddle point method: 



maxC(7r*7r*) ~ 2, maxC(7r*7r 

maxC(7r07rO) « 2 + 
According to the asymptotic expansion (|6^)-(|6 

2 



1 



G(t=0) 
G(t=0) 



(83) 



G{t=0) 
G{t=0) 



Git = 0) 



G{t = 0) 



G{t = 0) 



to/3 
(2TOr)3 
and 



for TOT ^ 1 



(84) 



G{t = 0) 



1 — (tttot)^ for mr ^ 1 



The corresponding plot of max C(mT) is demonstrated in Fig. 2. 
For pions this value is noticeable only at r wl fm and is 0.25 approximately, 
so max CttOttO (t = Ifm) w 2.25, max C7r-,r+ ("J" — Ifm) ~ 1.25. These values differ 
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from the quantum-mechanical resuhs pq ] where maxC^OTrO — 2, and C^^-T^^r = 1. 
This means that the effect has relativistic nature (decomposition of a field into the 
positive and negative frequency components and their interference in finite regions 
of homogeneity). 

Note that if the chemical potential \x tends to critical value /ic then even for 
large mr 3> 1 we have for maximal intercept (when the both quanta are very soft, 
Pi ~ P2 ~ 0) unusually large values: ~ 2 for tt+tt^ and ~ 3 for tt^tt" . For equally 
charged pions 7r+7r+ and tt^tt^ the intercept has the standard value ~2 for any mr 
and /i. 

In the typical experimental situation when mr 3> 1 and [i — Q the Bose-Einstein 
correlation functions have the standard structure (jl^) for all sorts of identical pions 
and longitudinal projection of the correlation function can be approximated by the 
expression 



exp 



C(p, <7i; qr = 0) « 1 + ^ ^ ^ ^ ^ "^'-^ 1 + exp [-r^AiglJ 

(l + T2A*gi)'' 

(85) 

where \l ^ \f^m^' ^^"^^ ^ "freeze out" temperature corresponding to proper 

time T when the particles leave the expanding matter. This result is obtained by 
the saddle-point method from ( ^0| ) and has the asymptotic form for longitudinal 
interferometry radius firstly obtained in Refs. for [imT ^ 1. 

The two- particle spectra for charge if +, ^-kaons are described by Eqs. (|80|), 
( ^ ) with substitution tt* because the complex representation for these 

fields can be replaced by the real one in the same manner as for charge pions. The 
same concerns of K^, K -pairs. Because of relatively large kaon mass the role of 
the addition term in the correlation function C{K^ , K^) is negligible. For the 
correlation functions of identical kaons asymptotic form in Eq.(p5|) can be used in 
all momentum region. 

The correlations in expanding photon gas are described by the formula (js^) with 
multiplier 1/2 at the second and third terms arising due to random polarization of 
photons. In this case the additional third term gives a good contribution for very 
soft photons producing maxC — 2 instead of =1.5 without the third term. 



7. Conclusions 

In this paper we give the theoretical analyses of spectra and correlations in inhomo- 
geneous weekly interacting boson gas. For the purpose the method of locally equilib- 
rium statistical operator used to calculate the averages such as {a'^{pi)...a{pj)...) has 
been developed. The problem was reduced to the system of the integro-differential 
and integral equations solved analytically for physically significant model of hydro- 
dynamic boost-invariant expansion. The main results are: 
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• the deviation of the particle phase-space density distribution from the Bose- 
Einstein one even in main approximation neglecting dissipative phenomena; 

• the appearance of the additional terms in the correlation functions of like 
and unlike (oppositely charged) particles as compared with the results of the 
nonrelativistic quantum mechanical approach. 

These effects are essential at small values rpo < 1 or/and at large enough 
chemical potentials /i ^ic. Under this condition the "effective" wave-length of 
the quanta, (po ~ m)^^ is larger than the length of homogeneity t in thcrmalizod 
medium and quanta begin to "feel" the all expanding matter. If = 0, it results 
in the additional number of soft quanta A\ic to interference of positive and negative 
frequency components of the relativistic quantum field in finite regions of homo- 
geneity. In special case of Bjorken boost-invariant picture this effect is described 
by the spectrum /m(po) of the so-called Milne's particles that appear at zero tem- 
perature in the hyperbolic space-time due to a mixing of the positive and negative 
frequency field components relatively to the Minkovski space. In the hydrodynamic 
picture the role of the Milne's Universe is played by the expanding thermalized 
matter "forming" this hyperbolic world by the isoterms. Here there are no effects 
at zero temperature: the spectrum is f{po) = fBE{po){'^fMiPo))- Note that Milne's 
spectrum /m cx |(aa)|^. The last value is responsible for the additional terms in the 
correlation functions. Therefore, the both effects have the common nature. They 
are connected with the space-time inhomogeneity of systems. 

The goal of the interferometry analysis in A + A collisions is to study the 
space-time evolution of the matter or, roughly speaking, to find proper time r 
of expansion. This is, actually, the average longitudinal length of homogeneity. 
Generally speaking, the "interferometry microscope" measures the size and shape 
of homogeneity regions in radiating sources. As known, the relative smallness of 
the effective emitting region is the basic condition for interferometry method to be 
applied experimentally. Under this circumstance the taking into account of the new 
effects for spectra and correlations of effectively soft bosons, {po — ^ t become 
to be important. 

The experimental consequences of the effect in the ultra-relativistic nucleus- 
nucleus collisions could also concern of the particles with small mass such as chiral 
quarks, gluons, photons, etc. The theory predicts the essential enhancement for 
a number of particles with small effective energy at an early stage of the matter 
expansion. This can lead to the increase in the number of photons and dileptons 
with small transverse momenta or small invariant masses if they are produced in 
the collisions of particles with small effective mass. 

Comments 

The paper is minor modified version of the unpublished preprint Yu.M.Sinyukov, 
ITP-93-8E, Kiev, 1993. Here was done firstly the interpretation of the HBT radii 
as the lengths of homogeneity in radiating systems. The modification takes into 
account some of the later results published in proceedings of the conferences: 

Yu.M.Sinyukov, Nucl. Phys. A566 (1994) 589c (QM 93); 
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Yu.M.Sinyukov. Spectra and correlations in small inhomogenious systems. 
In: Hot Hadronic Matter. Theory and Experiment, (J.Letessier, H.H.Gutbrod, 
J.Rafelski, eds.) p. 309, Plenum Publ., 1995. (NATO Workshop, Divonne-94); 

Yu.M.Sinyukov, S.V.Akkelin, A.Yu.Tolstykh, Nucl.Phys. A610 (1996) 278c 
(QM 96); 

Yu.M.Sinyukov, S.V.Akkelin, R.Lednicky, In Proc. of the 8th International 
Workshop on Maltiparticle Production in Matrahaza (T.Csorgo et al, eds), p. 66, 
World Scientific, 1998. 

and in the paper Yu.M.Sinyukov, B.Lorstad, Z.Phys.C61 (1994) 587. 
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Fig. 2. The aditional contribution to the correlation peak value for neutral and 
oppositely charged bosons at ^ = 0. 



